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HEAT  FLOW  IN  ATOMIZED  METAL  DROPLETS 

by 

★ 

C.  G.  Levi  and  R.  Mehrabian 


ABSTRACT 


The  solidification  of  spherical  droplets  with  a discrete  melt- 
ing temperature  is  analyzed  using  an  enthalpy  model.  Equations 
describing  the  cooling  of  the  initially  superheated  liquid  droplet 
and  a numerical  heat  flow  model  for  its  subsequent  solidification 
are  presented.  Important  parameters  like  times  for  initiation  and 
completion  of  solidification,  cooling  rates  and  interface  velocities 
in  aluminum,  iron  and  nickel  are  related  to  the  process  variables 
governing  the  rate  of  heat  extraction  from  the  droplets.  The  analy- 
sis is  performed  for  the  range  of  Biot  numbers  of  practical  interest 
where  Newtonian  cooling  models  are  not  considered  applicable,  0.01 
$ Bi  « 1.0,  and  the  results  are  presented  in  the  form  of  normalized 
or  dimensionless  quantities.  It  is  shown  that  the  average  cooling 
rate  in  the  liquid  prior  to  solidification  can  be  computed  with  the 
Newtonian  cooling  expressions.  However,  significant  temperature 
gradients  are  noted  at  the  droplet  surface  even  for  Biot  numbers  as 
low  as  0.01.  Reducing  the  droplet  diameter  reduces  the  times  nec- 
essary for  the  initiation  and  completion  of  solidification,  increases 

*C.  G.  Levi  and  R.  Mehrabian  are  Graduate  Student  in  the  Department 
of  Metallurgy  and  Mining  Engineering  and  Professor  in  the  Depart- 
ment of  Metallurgy  and  Mining  Engineering  and  the  Department  of 
Mechanical  and  Industrial  Engineering,  respectively.  University  of 
Illinois,  Urbana,  IL  61801. 
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\ the  interface  velocities  at  equivalent  fractions  solidified  and  de- 

creases  the  6,/R  ratio.  Although  smaller  droplet  diameters  promote 
higher  cooling  rates  in  the  liquid  at  the  beginning  and  in  the  solid 
i at  the  end  of  solidification,  the  effect  at  the  intermediate  stages 

is  more  complex  and  depends  on  the  initial  superheat,  the  Biot  num- 
ber  and  the  thermophysical  properties  of  the  material. 
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I.  INTRODUCTION 

During  solidification  of  metal  powders  heat  is  extracted  from 
the  droplets  by  both  convection  and  radiation  at  their  surface. 
However,  there  are  no  accurately  established  values  for  the  com- 
bined radiative  and  convective  heat  transfer  coefficient,  h,  and 
direct  measurement  of  the  cooling  rate  or  heat  flux  during  solidif- 
ication of  an  atomized  droplet  would  be  extremely  difficult,  if  not 

impossible.  In  gas  atomization  the  convective  heat  transfer  coef- 
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ficient  is  overriding  and  an  upper  limit  of  ^10  W/m  K can  be  esti- 
mated from  existing  expressions  for  h under  the  most  favorable  ex- 
perimental conditions  (1). 

Indirect  estimates  of  heat  transfer  coefficients  in  various 
atomization  processes  have  also  been  made  by  comparison  of  measured 
segregate  (dendrite  arm)  spacings  in  crystalline  alloy  powders  with 
predetermined  relationships  between  these  spacings  and  average  cool- 
ing rates  during  solidification.  The  heat  transfer  coefficients 
for  gas  atomization  deduced  in  this  way  are  of  the  same  order  of 
magnitude  as  those  estimated  from  mathematical  expressions  (1). 

In  general,  then,  a limitation  on  the  achievable  heat  transfer 
coefficient  at  a liquid  metal  droplet  - environment  interface  can 

be  translated  into  a limitation  on  the  important  dimensionless  var- 

★ 

iable,  Biot  number  , governing  the  rate  of  heat  extraction  from  the 

5 2 

droplet.  For  example,  a heat  transfer  coefficient  h ^ 10  W/m  .K 

_ ? 

translates  to  a range  of  Biot  numbers  of  10  < Bi  ^ 1.0  for  atom- 

ized droplets  of  liquid  aluminum  in  the  size  range  of  lym  to  lOOOyrn. 


* see  nomenclature 
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Previous  work  on  modelling  the  solidification  of  metal  drop- 
lets has  been  either  confined  to  the  Newtonian  cooling  regime  or 
to  a specific  range  of  heat  transfer  coefficients  (1-3).  In  this 
paper  general  relationships  are  developed  between  the  rate  of  heat 
extraction  and  the  important  solidification  parameters,  assumi ng 
the  heat  is  transferred  only  by  conduction  inside  a spherical  drop- 
let under  concentric  isotherms.  Since  Newtonian  cooling  expressions 
are  considered  applicable  for  Biot  numbers  below  .01  and  the  limit- 
ations on  h set  an  upper  limit  for  Bi  of  'v  1.0,  as  indicated  above, 
this  study  is  aimed  to  analyze  the  heat  flow  during  solidification 
of  superheated  metal  droplets  with  Biot  numbers  within  this  range. 
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II,  PROBLEM  STATEMENT  AND  SOLUTION  APPROACH 


We  consider  the  cooling  of  a superheated  pure  liquid  metal 
droplet  and  its  subsequent  solidification  due  to  radially  uniform 
heat  loss  by  radiation  and  convection  from  its  surface.  Surface 
conductance  - heat  transfer  coefficient  - and  the  temperature  of  the 
environment  are  assumed  to  remain  constant.  The  thermal  proper- 
ties of  the  solid  and  the  liquid  are  also  assumed  to  be  constant, 
but  different.  Heat  transfer  inside  the  droplet  is  by  conduction 
only. 

First,  a dimensionless  form  of  the  analytical  solution  to  the 
one-dimensional  heat  conduction  equation  in  spherical  coordinates 
is  used  to  determine  the  temperature  profile  in  the  droplet  when 
its  surface  reaches  the  melting  point.  This  profile  is  then  used 
as  the  initial  boundary  condition  for  the  subsequent  solidification 
problem.  The  concepts  of  the  enthalpy  model  recently  proposed 
by  Shamsundar  and  Sparrow  (4)  are  used  to  develop  the  numerical  sol- 
utions for  the  solidification  portion  of  the  problem.  In  this  model 
the  enthalpy  is  used  as  a dependent  variable  in  addition  to  the 
temperature,  and  the  finite  difference  equations  are  written  in  a 
fully  implicit  form.  The  important  variables  deling  cooling  and 
solidification  of  the  droplets  are  determined  and  expressed  in  their 
most  general  dimensionless  form  over  a range  of  Biot  numbers  of 
practical  interest  for  droplets  of  three  elements:  aluminum,  nickel 


'1 
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III.  MATHEMATICAL  DESCRIPTION 

For  mathematical  consistency  all  the  equations  in  the  follow- 
ing formulation  are  written  in  terms  of  the  dimensionless  Fourier 

* and  Biot  numbers  expressed  in  terms  of  the  properties  of  the  liquid. 

’ 

The  cooling  of  the  liquid  droplet  prior  to  solidification  was 
described  in  terms  of  the  analytical  infinite-series  solution  for 
the  radially  uniform  cooling  of  a sphere  with  a moderate  heat  trans- 
fer coefficient,  given  in  reference  (5). 

For  the  solidification  part  of  the  problem  the  spherical  drop- 
let is  divided  into  a fixed  number  of  concentric  shell  domains  of 
thickness  Ar  for  the  internal  nodes  (2<i<n),  and  Ar/2  for  the 
center  (i=l)  and  the  surface  (i  = n+1)  nodes.  An  energy  balance 
for  an  internal  volume  element  can  then  be  written  as: 
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where  a,.-,  a^  + and  v^  denote  the  inner  and  outer  surface  areas 
and  the  volume  of  the  spherical  nodal  shell  i,  respectively.  All 
the  variables  in  equation  (1)  and  those  given  below  are  defined  in 
the  Nomenclature. 

Equation  (1)  can  now  be  put  in  terms  of  dimensionless  nodal  en- 
thalpy, ^ , and  dimensionless  nodal  temperature,  6^ , and  becomes: 
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(3) 


For  the  center  node,  i = 1,  the  second  term  on  the  right 
hand  side  of  equation  (2)  becomes  zero. 

For  the  surface  node,  i = n + 1 , an  energy  balance  similar 
to  expression  (1)  is  written  which  includes  heat  loss  to  the  en- 
vironment. Appropriate  manipulation  of  this  equation  using  the 
definitions  of  dimensionless  nodal  enthalpy  and  temperature  yields 
the  following  corollary  expression  to  equation  (2): 


i 

|.j 

I 


n+1  3*  n+1  " ^n  + 1 jq  B 1 ^ 6 n + 1 + Ste^ 


In  general,  H..  in  equation  (1)  refers  to  the  specific  enthalpy 
of  the  volume  element  i and  assumes  different  forms  when  the  node 
is  in  the  liquid,  the  solid  or  contains  the  liquid-solid  interface, 
k in  equation  (3)  is  the  thermal  conductivity  of  the  liquid  or  the 
solid.  k+  = k~  for  all  nodes  except  that  containing  the  interface, 
i = £,  where  k+  = ks  and  k~  = kL. 

In  the  solid  region  the  dimensionless  nodal  enthalpy,  i|> . is 
always  negative  and  equal  to  the  dimensionless  nodal  temperature  0., 
whereas  in  the  superheated  liquid  region  \p.  is  larger  than  unity. 

A volume  element  containing  the  liquid-solid  interface  is  at  the 
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Equations  (2)  and  (4)  were  put  into  their  equivalent  finite 
difference  forms  using  backward  difference  formulas.  The  resulting 
system  of  simultaneous  equations  for  the  different  nodes  was  solved 
using  the  Gauss-Seidel  iterative  method  starting  with  the  tempera- 
ture distribution  obtained  from  the  analytical  solution  for  the 
pre-solidification  stage  (5).  The  convergence  criteria  used  for 
each  iteration  cycle  was: 


n + 1 

I 

i = 1 


t 

i . j + 1 


c 

i , j + 1 


) 


< 


( 5 ) 


where  the  superscripts  t and  c denote  the  trial  and  calculated  values, 
respectively. 

The  dimensionless  time  interval  AFo  was  adjusted  in  every  step 
j to  give  roughly  uniform  displacements  cf  the  liquid-solid  interface 
Trial  runs  were  performed  changing  the  number  of  nodes,  n+1,  and 
the  length  of  the  interface  displacements,  A<J>*,  (i.e.,  the  length 
of  the  time  steps  AFo).  It  was  observed  that  for  the  metals  inves- 
tigated, increasing  the  number  of  nodes  above  21  does  not  signifi- 
cantly affect  the  results.  On  the  other  hand,  some  parameters, 1 i ke 
Gl/R,  cooling  rates  in  the  liquid  and  interface  velocities  at 
the  beginning  (<J>*-+1)  and  at  the  end  ( <#>*-*-0 ) of  solidification  are 
particularly  sensitive  to  the  time  step  size,  especially  in  the  low 
range  of  Biot  numbers.  For  these  calculations  the  time  step  sizes 
were  appropriately  selected  until  convergence  was  reached.  Calcu- 
lations were  performed  on  a CDC  CYBER  175  compute*. 


IV.  RESULTS  AND  DISCUSSION 


The  analytical  equations  and  the  computer  heat  flow  model 
were  used  to  calculate  the  important  cooling  and  solidification 
parameters  (e.g.  cooling  rates,  interface  velocities,  tempera- 
ture gradients,  etc.)  as  a function  of  time  and  the  dimension- 
less variables  -mainly  Biot  Number-  governing  the  rate  of  heat 
extraction  from  the  droplets.  The  range  of  Biot  numbers  inves- 
tigated cover  the  values  of  practical  interest  in  liquid  metal 
atomization  where  Newtonian  cooling  models  are  not  considered 
applicable.  Calculations  were  carried  out  for  droplets  of  alum- 
inum, iron  and  nickel,  using  the  thermal  and  physical  properties 
listed  in  Table  I. 

In  the  first  part  of  this  section  the  results  obtained  from 
the  analytical  solution  for  the  cooling  of  the  liquid  are  described. 
Figure  1 shows  the  dependence  of  the  cooling  time,  t°,  (i.e.,  time 
for  the  surface  of  the  droplet  to  reach  its  melting  point)  on  the 
Biot  Number  and  initial  superheat.  The  time  calculated  from  the 
analytical  solution  has  been  normalized  for  purposes  of  represent- 
ation, dividing  it  by  the  time  obtained  from  the  Newtonian  cooling 
model,  which  is: 


Newt. 


In  (1  + ®o/Ste) 
~3Bi 


( 6 ) 


As  expected,  at  equivalent  initial  superheats  the  time  for 
non-Newtonian  cooling  is  shorter  because  not  all  the  sensible 
heat  in  the  droplet  is  removed  before  its  surface  reaches  the  melt- 
ing point.  Examples  of  calculated  cooling  times,  t°,  from  Figure  1 
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for  50  and  500pm  droplets  of  A1 , Fe  and  Ni  are  given  in  Table 
II  . Note  that  increasing  the  particle  diameter  by  one  order 
of  magnitude  (and  hence  the  Biot  number)  will  increase  the  cool- 
ing time  by  less  than  one  order  of  magnitude  as  indicated  by 
the  values  in  the  ordinate  of  Figure  1,  which  can  be  interpreted 
as  correction  factors  to  equation  ( 6). 

Normalized  temperature  distributions  in  liquid  droplets  when 
their  surface  reaches  the  melting  point,  for  different  Biot  Num- 
bers and  dimensionless  initial  temperatures  are  shown  in  Figure  2. 
The  largest  gradient  occurs  at  the  surface  and  is  given  by  the 
transformed  boundary  condition: 


) - Bi  * ( 


Ste< 


= Bi 
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( 7) 


It  should  be  noted  that  the  actual  temperature  gradient  at 
the  droplet  surface  is  independent  of  initial  superheat  and  drop- 
let size.  Furthermore,  calculations  show  that  the  effect  of  super- 
heat on  temperature  distributions  is  negligible  for  Bi  ~ .01. 
However,  the  results  also  indicate  that  even  for  small  Biot  num- 
bers, in  the  range  of  n,  0.01,  there  may  be  significant  temperature 
gradients  in  a metal  droplet.  For  example,  in  a 20pm  diameter  drop- 
let of  liquid  iron  where  Bi  ^ 0.01  (h  = 4 x 10^  W/m2K)  the  maxi- 
mum temperature  gradient  -at  the  droplet  surface-  when  t = t°  is 
n.  1.5  x 10®  K/m  and  a temperature  difference  of  7.5  K between 
the  surface  and  the  center  of  the  droplet  was  calculated  from 
the  analytical  solution. 
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Thus,  the  Newtonian  cooling  assumption  that  temperature  dif- 
ferences inside  a body  are  negligible  for  Bi  ~ 0.01  may  not  always 
be  justified.  Other  examples  of  gradients  at  the  surface  of  various 
droplets  are  given  in  Table  II. 

An  important  variable  affecting  supercooling  prior  to  nuclea- 
tion  is  the  cooling  rate  in  a liquid  droplet.  The  instantaneous 
average  cooling  rate  ft  estimated  from  the  series-type  solution  was 
found  to  be  identical  to  that  obtained  from  the  Newtonian  cooling 
approximation.  When  the  surface  of  the  droplet  reaches  its  melting 
point  (t  = t°)  the  cooling  rate  is: 


8.°  = 
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Examples  of  the  calculated  cooling  rates  for  50pm  and  500pm  drop- 
lets of  A1 , Fe  and  Ni  assuming  a high  heat  transfer  coefficient  of 
5 x 104  W/m^  K are  given  in  Table  II  . 

The  next  part  of  this  section  contains  the  results  obtained 
from  the  computer  simulation  of  heat  flow  during  solidification 
using  the  temperature  distributions  shown  in  Figure  2 as  initial 
conditions. 

The  net  solidification  time  for  droplets  of  A1 , Fe  and  Ni 
as  a function  of  Biot  number  and  dimensionless  initial  superheat 
is  presented  in  Figure  3.  The  time  resulting  from  the  numerical 
solution  of  the  enthalpy  model  has  been  normalized  dividing  it 
by  the  time  calculated  from  the  following  expression  for  solidif- 
ication in  a Newtonian  regime: 
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As  anticipated,  increasing  initial  superheat  prolongs  solid- 
ification, and  the  effect  is  larger  as  the  Biot  Number  increases; 
the  sensible  heat  retained  in  the  liquid  portion  of  the  droplet 
increases  resulting  in  longer  times  from  initiation  to  completion 
of  solidification.  For  example,  500ym  droplets  of  Fe  solidifying 

with  an  h = 5 x 10^  W/m^K  (Bi  = 0.31)  will  require  4.8  x 10"^s  to 
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complete  solidification  with  no  initial  superheat  and  5.3  x 10  s 
with  a dimensionless  initial  superheat  of  0.1  1 50K ) , see  Table  II 

The  effect  of  superheat  on  net  solidification  times  diminishes  with 
decreasing  Biot  Numbers. 

Figure  4 shows  dimensionless  liquid-solid  interface  velocities 
as  a function  of  Biot  Number  for  various  positions  of  the  interface 
in  A1  droplets  solidifying  with  concentric  isotherms.  The  data 
show  that  dimensionless  interface  velocity  increases  with  Biot  Num- 
ber and  the  progress  of  solidification. 

The  dimensionless  interface  velocity  for  Newtonian  solidif- 
ication is: 


Newt. 


Bi 


(10) 


For  a given  interface  position,  <J>*,  equation  (10)  would 
give  a straight  line  with  a slope  of  unity  In  Figure  4.  The  ratio 
of  calculated  non-Newtonian  to  Newtonian  interface  velocities  is 


■ 4. 


close  to  unity  in  the  beginning  of  solidification  for  the  case  of 
no  superheat  regardless  of  Biot  number.  On  the  other  hand,  this 
ratio  decreases  with  increasing  fraction  solid  and  increasing  Biot 
number  at  a given  fraction  solid.  In  general,  increasing  initial 
superheat  in  the  non-Newtonian  regime  results  in  a corresponding 
decrease  in  solid-liquid  interface  velocites.  However,  this  effect 
is  only  noticeable  in  the  initial  stages  of  solidification,  and 
decays  rapidly  with  increasing  superheat  and  decreasing  Biot  num- 
ber. 

Examples  of  calculated  interface  velocities  are  given  in  Table 
II.  Note  that  during  non-Newtonian  solidification  increasing  the 
particle  size  by  one  order  of  magnitude  (without  changing  h)  re- 
duces the  actual  interface  velocity  at  a particular  fraction  solid; 
On  the  other  hand,  increasing  the  heat  transfer  coefficient  h 
(without  changing  r )will  increase  R at  any  location. 

The  results  of  normalized  ratio  of  temperature  gradient  in  the 
liquid  at  the  liquid-solid  interface,  GL,  to  interface  velocity  R 
for  the  three  metals  and  a dimensionless  initial  superheat  8Q/Ste 
of  0.3  are  presented  in  Figure  5.  GL  at  the  beginning  of  solidifi- 
cation is  given  by  the  boundary  condition  at  the  surface: 


G 


L 


$*»1 


-r-  Bi  Ste  • 
r _ 


(11) 


while  the  initial  interface  velocity  is  close  to  that  calculated 


_ . 
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from  equation  (10),  hence,  the  initial  GL/R  ratio  should  asympto- 
tically approach  a constant  value  with  decreasing  Biot  number: 


SL 

C7 


when  Bi  ■+■  0 


(12) 


Examples  of  GL/R  ratios  calculated  from  Figure  5 are  listed  in 
Table  II.  As  expected,  this  parameter  decreases  rapidly  with  the 
progress  of  solidification,  and  with  Biot  number  for  a given  posi- 
tion of  the  liquid-solid  interface.  Furthermore,  reduction  of  the 
initial  superheat  (not  shown  in  Figure  5)  results  in  a corresponding 
decrease  in  GL/R,  the  effect  being  more  pronounced  when  0o/Ste  is 
in  the  range  of  0 to  0.1. 

Normalized  values  of  the  calculated  average  cooling  rates  in 
the  solid  and  the  liquid  during  solidification  of  iron  droplets  are 
shown  in  Figures  6 and  7,  respecti vely.  Note  that  the  part  of  the 
normalizing  factor  inside  the  brackets  is  the  Newtonian  net  solid- 
ification time  obtained  from  equation  ( 9),  which  is  directly  pro- 
portional to  rQ  and  inversely  proportional  to  h.  Therefore,  in- 
creasing the  Biot  number  by  increasing  the  droplet  size,  rQ,  will 
increase  the  actual  cooling  rate  only  if  the  slope  of  the  curve  at 
that  point  is  larger  than  unity.  On  the  other  hand,  increasing  the 
heat  transfer  coefficient  h will  always  increase  the  cooling  rate 

The  average  cooling  rate  in  the  solid,  ?2S,  Figure  6,  increases 
with  the  progress  of  solidification  for  Bi  < 0.9.  However,  the 
results  also  indicate  that  for  higher  Biot  numbers  the  cooling  rate 
has  a maximum  at  some  intermediate  stage  during  solidification, 

(((>*  f 0).  Finally,  decreasing  the  initial  superheat  will  always 
increase  the  achievable  cooling  rates  in  the  solid,  but  the  effect 
is  only  felt  at  the  initial  stages  of  solidification. 


- -*  V.  . 
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The  average  cooling  rate  in  the  liquid,  ftL,  for  iron.  Figure 
7 , decreases  with  the  progress  of  solidification;  the  effect  being 
more  pronounced  at  small  Biot  numbers. 

Examples  of  average  cooling  rates  JT<.  and  calculated  from 
Figures  6 and  7 and  similar  plots  obtained  for  aluminum  and  nickel 
(showing  identical  trends)  are  given  in  Table  II  . 
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SUMMARY 


The  infinite  series  analytical  solution  to  the  cooling  of  an 
initially  superheated  metal  droplet  under  concentric  isotherms  was 
coupled  to  a numerical  solution  for  the  heat  flow  problem  during 
the  subsequent  solidification,  written  in  terms  of  both  enthalpy 
and  temperature  as  dependent  variables.  The  results  obtained  from 
this  analysis  as  compared  to  those  predicted  by  the  Newtonian  cool- 
ing model  for  aluminum,  iron  and  nickel  are  summarized  below: 

1)  Both  models  predict  that  the  average  cooling  rate  in  the  liquid 
prior  to  solidification,  is  directly  proportional  to  the 
heat  transfer  coefficient,  h,  inversely  proportional  to  the 
droplet  radius  r , and  independent  of  initial  superheat. 

2)  The  cooling  time  (for  the  surface  to  reach  the  melting  point), 
t°,  estimated  from  the  analytical  solution  increases  with  in- 
creasing superheat  towards  the  Newtonian  model  prediction  (which 
is  a logarithmic  relationship).  Reducing  rQ  has  a less  pro- 
nounced effect  on  t°  than  the  direct  proportionality  indicated 

by  the  Newtonian  model.  On  the  other  hand,  increasing  h has  a 
more  pronounced  effect  on  decreasing  t°  than  the  inverse  propor- 
tionality predicted  by  the  Newtonian  expression. 

3)  Substantial  temperature  gradients  at  the  surface  of  the  droplet 
when  t = t°  are  calculated  from  the  analytical  solution  for  Biot 
numbers  as  low  as  0.01.  Although  these  gradients  are  independent 
of  superheat,  the  latter  markedly  influences  the  internal  temp- 
erature profiles  for  Bi  Z 0.01. 

4)  Reducing  rQ  diminishes  the  net  solidification  time,  t^,  in  a 
more  pronounced  way  than  the  direct  proportionality  predicted 
by  the  Newtonian  model.  Furthermore.it  increases  the  interface 


velocity  R at  equivalent  fractions  solidified  towards  the 
values  predicted  by  the  Newtonian  model  (which  are  independ- 
ent of  r ).  Reducing  rQ  also  decreases  the  G^/R  ratios  and 
increases  the  achievable  cooling  rates  in  the  liquid  at  the 
beginning  of  solidification,  4>*-*-l  , and  in  the  solid  at  the 
end  of  it,  4>*-*0.  The  effect  of  rQ  on  iTL  and  iTs  at  intermed- 
iate stages  during  solidification  is  more  complex  and  depends 
on  the  metal,  the  corresponding  Biot  numbers  and  the  initial 
superheat,  as  well  as  the  fraction  solidified. 

5)  Increasing  h reduces  t and  increases  R in  a less  pronounced 
manner  than  that  predicted  by  the  Newtonian  model.  It  always 
increases  GL/R  as  well  as  the  cooling  rates  and  at  any 
stage  during  solidification.  The  progress  of  solidification 
always  increases  the  interface  velocity  and  decreases  GL/R  as 
well  as  cooling  rates  in  the  liquid.  However,  the  cooling  rate 
in  the  solid  increases  monotonical ly  with  time  at  the  smaller 
Biot  numbers  but  apparently  shows  a maximum  before  the  end  of 
solidification  above  a certain  Biot  number,  which  for  iron  is 
~0.9  and  for  aluminum  and  nickel  is  above  1.0. 

6)  Increasing  initial  superheat  increases  tf,  GL/R  and  ?TL,  but  de- 
creases R and  The  effects  of  superheat  are  mostly  felt  in 

the  initial  stages  of  solidification  and  decay  rapidly  with  de- 
creasing Biot  number. 


NOMENCLATURE 


Dimensioned  variables 


a 


i 


C 

G 

h 

H 


r 


r 

o 

R 

t 

t° 


Tg=300°K 


M 


vi 


a=  k/pC 
ft 


P 


area  of  the  inner  (-)  or  outer  (+)  boundary 
surface  of  the  shell  element  i 
specific  heat 

gradient  in  the  liquid  at  the  liquid-solid  interface 
heat  transfer  coefficient  at  the  droplet  surface 
specific  enthalpy 

specific  enthalpy  of  the  solid  at  TM  (reference  value), 
heat  of  fusion 
thermal  conductivity 

position  along  the  radial  coordi nate,  0 < r < rQ . 
radius  of  the  spherical  droplet 
interface  velocity 
time 

cooling  time 
net  solidification  time 
temperature 

temperature  of  the  cooling  medium 
initial  temperature  of  the  droplet, 
melting  temperature 
volume  of  the  shell  element  i 
thermal  diffusivity 
instantaneous  average  cooling  rate 
density,  average  of  liquid  and  solid  at 
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D i me n s i o n 1 ess  variables 


Bi  = 

hr0/kL 

Biot  number 

II 

O 

Li. 

aLt/ro 

Fourier  number  (or  time) 

AFo  = 

“L4t/ro 

time  increment 

Ste  = 

cs  (tm-tg)/4hsl 

Stefan  number 

$ = r/r. 


r*/r 


A<}>  = Ar/r, 


♦i  • *i 


0 = C$(T  - Tm)/AHsl 


9o  = CS^To"TM^/AHSL 


eQ/Ste 


H>  = (H  - Hs)/AHsl 


posi tion 


fractional  position  of  the  liquid- 
solid  interface 


nodal  spacing 


radius  of  the  surfaces  at  r.  - Ar/2 
and  ri  + Ar/2,  respecti vely ;for  the 

center  node  <J>7  = 0 and  for  the  surface 

node  <t>*+1  = 1 . 


temperature 
initial  temperature 
superheat 
enthalpy 


- 

i 


»-  •--  - -- 


_ ...  _ /.  . ...  , 


f 

i 

► 


f = denotes  completion  of  solidification 

i - nodal  index,  1 £ i s n + 1,  i = i for  r = 0,  i = n + 1 for  r=r 
j = time  step  index,  j > 1 

l = index  of  volume  element  containing  the  liquid-solid 
interface 

n = number  of  space  intervals  between  nodes 


G = denotes  gas  environment 
L = liquid 
M = denotes  melting 
S = solid 


Superscr i pts 


+ in  the  positive  r direction 

in  the  negative  r direction 
* liquid-solid  interface 

° denotes  the  moment  when  the  surface  reaches  the  melting 
temperature  and  solidification  starts 
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Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Normalized  time  for  the  surface  of  a liquid  droplet  to  cool 
from  an  initial  temperature  of  Tq  to  its  melting  temperature 
T^  by  convection  to  an  environment  with  constant  temperature 
Tg,  as  a function  of  the  Biot  number  and  dimensionless  initial 
superheat  0o/Ste. 

Normalized  temperature  distributions  in  a liquid  droplet 
when  its  surface  reaches  the  melting  temperature  TM,  for 
different  Biot  numbers  and  dimensionless  initial  superheats 

Vste  = (To  " tm)/(tm  - V*  I 

Normalized  net  solidification  time  for  liquid  droplets  of 

aluminum,  iron  and  nickel,  as  a function  of  Biot  number  and 

dimensionless  initial  superheat  0o/Ste.  The  calculated  curve 

for  nickel  with  no  superheat  closely  follows  that  of  aluminum 

for  0.  = 0.3Ste. 

0 

Dimensionless  liquid-solid  interface  velocity  during  solidifi- 
cation of  aluminum  droplets,  as  a function  of  Biot  number  and 
fractional  position  of  the  interface  4>*. 

Normalized  ratio  of  the  temperature  gradient  in  the  liquid 
at  the  liquid-solid  interface,  G^,  to  the  interface  velocity, 

R,  during  solidification  of  A1 , Fe  and  Ni  droplets,  as  a 


function  of  Biot  number  and  fractional  position  of  the  inter- 
face. 0q  and  Ste  denote  the  dimensionless  initial  temperature 
and  Stefan  number,  respectively. 


r 
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Figure  6 Normalized  instantaneous  average  cooling  rate  in  the  solid 

during  solidification  of  iron  droplets,  as  a function  of 
Biot  number  and  fractional  position  of  the  solid-liquid  inter- 
face 4>*.  0 and  Ste  denote  the  dimensionless  initial  temper- 

ature and  Stefan  number,  respectively. 

Figure  7 Normalized  instantaneous  average  cooling  rate  in  the  liquid 

during  solidification  of  iron  droplets,  as  a function  of  Biot 
number  and  fractional  position  of  the  solid-liquid  interface 
<J>*=r*/rQ.  60  and  Ste  denote  the  dimensionless  initial  temp- 

erature and  Stefan  number,  respectively. 
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TABLE  I 

PROPERTIES  OF  METALS  USED  IN  THE  CALCULATIONS* 


Properties 

A1 

Fe 

Ni 

Melting  Temperature  TM  (K) 

933 

1810 

1726 

Heat  of  Fusion  x 10  5 

3.95 

2.72 

3.02 

__  O O 

Average  Density  p (Kg/m  ) x 10 

2.46 

7.19 

8.07 

Specific  Heat  ( J / Kg  K)  x 10-3 

Liquid  C^ 

1.09 

0.787* 

0.734 

Solid  Cs 

1.19 

0. 764+ 

0. 61 9+ 

Thermal  Conductivity  (W/mK) 

Liquid  kL 

90.8 

40.3 

39. 3+ 

Solid 

210.8 

34.6 

87.0  + 

Thermal  Diffusivity  a=k/pC  (m^/s)xl0+^ 

Liquid  aL 

3.40 

0.712 

0.664 

Solid 

7.18 

0.630 

1.74 

Stefan  Number  (for  Tg  = 300K) 

1.91 

4.24 

2.93 

Actual  Superheat,  Tq-Tm  f°r 

To'TM  °‘4 

253 

604 

570 

Dimensionless  0Q/Ste  = j— _-y - = q 3 

MG 

190 

453 

428 

0.2 

127 

302 

285 

0.1 

63 

151 

143 

*Thermophysical  properties  at  T^ 
^Estimated  by  extrapolation 


[3Bi/ln(60/Ste  + 1)] -(aLt°/r0 


BIOT  NUMBER  (Bi  = hrQ/kL) 


Figure  1.  Normalized  time  for  the  surface  of  a liquid 
droplet  to  cool  from  an  initial  temperature 
of  Tq  to  its  melting  temperature  TM  by  con- 
vection to  an  environment  with  constant  temp 
erature  Tq,  as  a function  of  the  Biot  number 
and  dimensionless  initial  superheat  e0/Ste. 
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Figure  2.  Normalized  temperature  distributions  in  a 
liquid  droplet  when  its  surface  reaches 
the  melting  temperature  T^,  for  different 
Biot  numbers  and  dimensionless  initial 
superheats  eQ/Ste  = (to-Tm^ ^TM'TG^  * 


NORMALIZED  SOLIDIFICATION  TIME  (3-Bi -Ste) (CL/C$)- (aLtf/rQ2) 


NET  SOLIDIFICATION  TIME 


cs(W 


Figure  3. 


BIOT  NUMBER  (Bi  = hrQ/kL) 


Normalized  net  solidification  time  for  liquid 
droplets  of  aluminum,  iron  and  nickel,  as  a 
function  of  Biot  number  and  dimensionless  ini 
tial  superheat  e0/Ste.  The  calculated  curve 
for  nickel  with  no  superheat  closely  follows 


that  of  aluminum  for 


0. 3Ste . 


i 

Figure  5.  Normalized  ratio  of  the  temperature  gradient 
in  the  liquid  at  the  liquid-solid  interface, 

Gi  , to  the  interface  velocity,  R,  during 
solidification  of  A1 , Fe  and  Ni  droplets,  as 
a function  of  Biot  number  and  fractional  posi- 
tion of  the  interface  6<j  and  Ste  denote  the 
. dimensionless  initial  temperature  and  Stefan 

number,  respectively. 
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NORMALIZED  AVERAGE  COOLING  RATE  IN  THEJSOLID 
[(r*/a,  )(CQ/C.  )/(3BiSte)](C./AH<-1  )fl- 


I 


Figure  6.  Normalized  instantaneous  average  cooling 
rate  in  the  solid  during  solidification 
of  iron  droplets,  as  a function  of  Biot 
number  and  fractional  position  of  the 
solid-liquid  interface  4>*.  0Q  and  Ste 
denote  the  dimensionless  initial  temper- 
ature and  Stefan  number,  respectively. 
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NORMALIZED  AVERAGE  COOLING  RATE  IN  THE_L IQUI D 
[(r‘/a. )(C«./C,  )/ ( 3Bi Ste ) ] (C,  /AH-.  )ft, 


10"2  10_1  10° 
BIOT  NUMBER  (Bi  = hrQ/kL^ 


Figure  7.  Normalized  instantaneous  average  cooling 
rate  in  the  liquid  during  solidification 
of  iron  droplets,  as  a function  of  Biot 
number  and  fractional  position  of  the 
solid-liquid  interface  <\>*=r*/r0.  e0  and 
Ste  denote  the  dimensionless  initial  temp- 
erature and  Stefan  number,  respectively. 


RAPID  MELTING  AND  SOLIDIFICATION  OF  A SURFACE 
DUE  TO  A STATIONARY  HEAT  FLUX 


by 

S.  C.  Hsu,  S.  Kou  and  R.  Mehrabian* 


ABSTRACT 

A general  two-dimensional  computer  heat  flow  model  is  devel- 
oped in  an  oblate  spheroidal  coordinate  system  for  rapid  melting  and 
subsequent  solidification  of  the  surface  of  a semi -inf i ni te  solid 
subjected  to  a high  intensity  heat  flux  over  a circular  region  on 
its  bounding  surface.  Generalized  numerical  solutions  are  presented 
for  an  aluminum  substrate  subjected  to  both  uniform  and  Gaussian 
heat  flux  distributions.  Temperature  distributions,  melt  depth 
and  geometry,  and  melting  and  solidification  interface  velocities 
are  calculated  as  a function  of  applied  heat  flux,  radius  of  the  cir- 
cular region,  and  time.  It  is  shown  that  the  important  melting 
and  solidification  parameters  are  a function  of  the  product  of  the 
absorbed  heat  flux,  £,  and  the  radius  of  the  circular  region, 

General  trends  established  show  that  melt  depth  perpendicular  to 
the  surface  is  inversely  proportional  to  the  absorbed  heat  flux 
for  a given  temperature  at  the  center  of  the  circular  region.  Dimen- 
sionless temperature  distributions  and  the  ratio  of  liquid-solid 
interface  velocity  to  absorbed  heat  flux,  R/q , as  a function  of 
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dimensionless  melt  depth  remain  the  same  if  the  product  £a  is  kept 
constant,  while  £ and  £ are  varied.  For  a given  total  power, 
absorbed  melting  and  solidification  parameters  are  compared  for 
uniform  and  Gaussian  heat  flux  distributions.  For  a given  tempera- 
ture at  the  center  of  the  circular  region  both  melt  depth  and  width 
are  smaller  for  the  Gaussian  distribution  while  temperature  grad- 
ients and  interface  velocities  are  larger. 


I.  INTRODUCTION 


In  a recent  paper  rapid  melting  and  solidification  of  a met- 
allic surface  layer  subjected  to  high  intensity  heat  fluxes,  such 
as  those  achieved  via  continuous  CO^  laser  radiation,  was  dis- 
cussed (1).  The  effect  of  absorbed  heat  flux  on  the  important 
surface  layer  melting  and  subsequent  solidification  variables  of 
three  substrate  materials;  aluminum,  iron  and  nickel  were  deter- 
mined using  a simplified  one-dimensional  computer  heat  flow  model. 

A second  limitation  imposed  on  the  previous  problem  statement  was 
that  only  uniform  heat  fluxes  were  considered.  In  general,  both 
of  these  restrictions  are  violated  in  practice.  Exception  is  noted 
for  cases  where  the  spot  diameter  of  the  incident  radiation  is 
large  relative  to  the  heat  affected  zone  in  the  work  piece  and 
the  one-dimensional  heat  flow  model  for  the  case  of  uniform  heat 
flux  is  applicable. 

This  paper  addresses  the  two-dimensional  (axi-synmetrlc) 
heat  flow  problem  of  the  melting  and  subsequent  solidification  of 
a surface  layer  subjected  to  a stationary  high  intensity  heat  flux 
over  a circular  region  on  its  boundary  surface.  A general  method 
for  the  solution  of  this  type  of  a problem  is  developed  which  per- 
mits treatment  of  both  uniform  and  non-uniform  heat  flux  distribu- 
tors in  the  circular  region.  Results  of  numerical  computations 
carried  out  on  a Control  Data  CYBER  175  computer  are  presented  for 
a pure  aluminum  substrate  subjected  to  both  uniform  and  Gaussian 
heat  flux  distributors. 

The  two  main  assumptions  of  this  paper  are:  (a)  thermal  proper- 
ties of  the  liquid  and  the  solid  are  independent  of  temperature,  but 
different  from  one  another,  and  (b)  convection  in  the  liquid  region 
is  not  taken  into  consideration. 


II.  PROBLEM  STATEMENT  AND  SOLUTION  APPROACH 


Consider  the  rapid  melting  and  subsequent  solidification 
of  the  surface  layer  of  a semi  - inf i ni te  solid,  initially  at  room 
temperature,  subjected  to  a high  intensity  heat  flux  over  a cir- 
cular region  on  its  bounding  surface.  The  surface  outside  the 
heated  region  is  considered  adiabatic.  The  heat  conduction  equa- 
tion for  the  ax i -symmetri c system  is  put  in  oblate  spheroidal 
coordinates  which  is  a "more  natural"  coordinate  system  for  this 
problem  geometry,  i.e.,  is  more  compatible  with  the  bounding  surface. 
This  results  in  reduced  storage  requirements  as  well  as  computa- 
tional time.  The  generalized  expressions  derived  by  Yovanovich 
et  al  ( 2 ) for  the  determination  of  the  coefficients  in  the  fin- 
ite difference  equations  governing  the  heat  transfer  within  the 
discretized  spatial  domain  are  used.  However,  in  contrast  to 
previous  treatments  of  melting  and  solidification  problems  the 
finite  difference  equations  are  rewritten  in  terms  of  the  enthalpy 
model  recently  proposed  by  Shamsundar  and  Sparrow  (4)  for  multi- 
dimensional phase  change  problems.  In  this  model  the  enthalpy  is 
used  as  a dependent  variable  in  addition  to  the  temperature.  In 
this  way  a single  energy  conservation  equation  is  used  for  each 
discritized  spatial  domain  which  by  appropriate  substitution  of 
physical  constants  is  applicable  regardless  of  whether  the  domain 
is  in  the  liquid  state,  the  solid  state  or  contains  the  liquid- 
solid  interface.  The  interface  shape,  its  velocity  and  other  im- 
portant melting  and  solidification  parameters  are  obtained  by  per- 
forming auxiliary  computations.  Finally,  the  computational  tech- 
nique utilized  is  a fully  implicit  one. 
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III.  MATHEMATICAL  DESCRIPTION 


The  axi -symmetric  oblate  spheroidal  coordinate  system  used 
in  this  study  is  shown  in  Figure  1.  The  relationships  between 

the  circular  cylinder  and  the  oblate  spheroidal  coordinates  are: 
r = a cosh  n sin  C ^1) 

z = a sinh  n cos  £ 

where  a^  is  the  radius  of  the  circular  region  subjected  to  the 

heat  flux. 

The  control  volume  is  centered  about  an  arbitrary  point  (i,j) 
in  space  in  Figure  1,  where  i and  j indicate  the  finite  discriti- 
zation  of  space  in  the  n and  £ directions,  respectively.  The 
change  in  respective  coordinate  values  between  successive  nodes 
in  the  n and  £ directions  are  An  and  A£,  respectively. 

The  general  formulation  of  the  finite  difference  represent- 
ation of  the  heat  conduction  equation  in  orthogonal  curvilinear 
coordinates  and  its  specific  form  in  the  discritized  domain  of 
Figure  1 centered  about  node  (i,j)  is  given  by  Yovanovich  et  al 
(2,3).  For  the  specific  problem  at  hand  this  equation  can  be  put 
in  the  following  backward  difference  form: 

mm  mm  mm 

-C1  <TM  - * C2  <Vl.J  - T1.J>  - C3  <Ti,j  - 


+ C.  (T, 


i , j + 1 


m 

T,-  ,)  + 


2irk  r j sA£An 

5 I » J 


cosh2ni  - s i n 2 £ . m m_i 

m 1 <Ti.j  - Ti  .j  > 

Where  the  coefficients  are: 

, tanh  Hi 


, cot  £ 
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k At  a At 

AFo  = — * — — 

pCpsa2  a2  (3)  continued 

r.  . = a cosh  n,  sin  £. 

1 * J I J 


and  PAV  is  the  total  source  strength  of  the  arbitrary  control  vol- 
ume. While  the  subscripts  i and  j refer  to  spatial  location,  the 
super  subscript  m refers  to  time  level.  All  the  other  terms  in 
the  equations  are  defined  in  the  Nomenclature. 

The  finite  difference  representation  of  the  heat  conduction 
equation,  expression  (2),  can  now  be  put  in  terms  of  dimensionless 
nodal  enthalpy,  ip , and  dimensionless  nodal  temperature,  6.  These 
two  dependent  variables  are  defined  as: 


* ' sW/ 


(H  - H ) 


0 = 


JL  ^£S. 


*Hs£ 

<T  - V 


dV  = 


H - H 


AH 


S£ 


AH 


Si 


(4) 


In  general  H in  equation  (4)  refers  to  the  specific  enthalpy 

of  the  discritized  space  volume  and  assumes  different  forms  when 

the  node  is  in  the  liquid,  solid  or  in  the  liquid-solid  region. 

★ 

H$  is  the  specific  enthalpy  of  the  solid  at  its  melting  point  and  k 
is  thermal  conductivity  which  assumes  conductivity  of  the  liquid 
or  the  solid.  In  the  formulation  presented  here  the  following 
assumptions  are  made;  the  density  of  the  solid  and  the  liquid  are 
equal,  the  conductivities  and  heat  capacities  of  the  solid  and  li- 
quid are  independent  of  temperature,  but  different  from  each  other. 


In  the  solid  region  the  dimensionless  nodal  enthalpy,  \p,  is 
negative  and  is  equal  to  the  dimensionless  nodal  temperature,  6: 


ip  = 6 = 


^s_(T  - V 


AH 


< 0 


s Jt 


In  the  superheated  liquid  region: 
Cp,(T  - V 

* ■ ' ♦ — - > '•» 


s % 


0 = 


k S CP£ 


U>  - 1) 


For  an  element  containing  the  liquid-solid 
ing  point  of  the  metal,  ip  is  defined  by: 


(5) 


(6) 


interface  is  at  the  melt 


0 < ip  $ 1.0  and  0=0  (7) 

The  value  of  \p  is  equal  to  the  weight  fraction  of  the  element 
which  is  in  the  liquid  state,  f^. 

The  appropriate  forms  of  equation  (4)  for  an  element  in  the 
liquid  or  solid  phase  are  substituted  into  equation  (2).  After 
some  manipulation  and  multiplication  by  different  factors  the  fol- 
lowing common  equation  is  obtained  for  both  cases. 


cosh2  n.  - sin2  £.  m m 

AFo  1 *i,j  + (C1  + C2  + C3  + C4}  ®i,J  = 


cosh2  n-  - sin2  m-1  m 

- J 


m 


AFo 


*1.j  + (C1  ®i-l,j  + C2  ®1+l,j  + C3  ®1,j-l 
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(8) 
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IV.  BOUNDARY  CONDITIONS  AND  SOLUTION  OF  THE 
FINITE  DIFFERENCE  EQUATIONS 


The  Initial  and  boundary  conditions  are  derived  based  on  the 
following  assumptions. 

The  metal  is  initially  at  room  temperature , T = T . 


The  absorbed  heat  flux  in  the  circular  region  of  the 
surface  is  in  general  q = q(r)  and  the  problem  is 
solved  for  the  two  special  cases  of  constant  and 
Gaussian  distribution  of  the  absorbed  heat  flux. 


Equation  ( 8 ) is  then  subject  to  the  following  boundary  con- 
d i tions : 


(i) 

t = 0, 

o^n<®»  o x<  ( s<  ' 

i .j 

» CPS  <tm  - T 

i.j--  ahsA 

(ii) 

t > 0, 

n = 0,  0 < C < TT/2 

(a)  q 

= Q/ua2  = constant 

(b)  q 

_ -2  sin2C 

= q0  e 

(iii) 

t > 0, 

£ = tt/2,  n > 0 

o 

4* 

II 

O 

(iv) 

t > 0, 

n oo,  0 4:  c ^ Tf/2 

(9  ) 
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(11) 
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(v)  t > 0,  n > 0,  5 - 0 


C3  = 0 


(13) 
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Boundary  condition  (iii)  assumes  that  convective  and  radiative 
heat  loss  from  the  surface  are  negligible.  This  is  in  line  with 
our  earlier  calculations  (1). 


- -rl  t tii4L  — - I.  v 
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Finally,  for  the  boundary  conditions  (ii)  to  (v)  the  value  of 
PAV  in  equation  (8)  is  zero  except  at  the  circular  region  on  the 
surface,  boundary  condition  (ii),  when  it  is  specified  by  the  fol- 


lowing equations: 


C.  + AC/2 

PAV  = ira2  f J q (€)  sin  (2£)  d£  (14) 

«J  - M/2 


For  case  (a)  q = constant 


PAV  = 2ira2q  sin?.  cos£.  sinA£ 

J J 

For  case  (b)  q = qQe  "2sin  5 
PiV  - "‘ip  [ e'2  si"!  (c  - 


2 sin2  (C  + ^)-j 


The  system  of  algebraic  equations,  equation  (8)  for  the  dif- 


ferent elements  in  the  oblate  spheroidal  coordinate  were  solved 


using  the  Gauss-Seidel  iterative  method.  The  computer  logic  pre- 
sented below  closely  follows  that  described  by  Sparrow  et . al.  (4  ) 
However,  the  family  of  equations  and  specific  relationships  between 

dimensionless  nodal  enthalpy  and  temperature  developed  here  are 
different  in  problem  geometry,  coordinate  system,  boundary  con- 
ditions and  the  fact  that  superheat  in  the  liquid  is  taken  into 


consideration . 


Once  the  solution  las  been  obtained  for  a time  level  m - 1, 
the  calculation  for  tha  following  time  level  m is  started  by 


initially  setting  all  the  . values  as  equal  to  the  previous 
time  level.  Then  using  the  equations  (5)  to  (7)  the  left-hand 
side,  L.H.S.,  of  equation  (8)  is  written  in  terms  of  dimensionless 
nodal  enthalpy  as  follows;  If 'J'fj  < 0 *he  element  is  in  the  solid. 


♦ i*  • * 07  *,and  both  the  left  and  right-hand  sides  of  equation 

(8)  are  less  than  zero.  On  the  other  hand,  if  0 < Vj1  < < 1.0, 

1 *J 

the  element  contains  the  liquid-solid  interface  and  the  following 
relationship  holds  for  the  R.H.S.  of  equation  (8). 


0 < R.H.S.  < 


cosh2n-  - sin2£, 


Finally,  if  ^ . > 1.0  the  element  is  in  the  superheated  liquid 
• » J 

region  and  substitution  of  equation  (6)  into  equation  (8)  yields: 

cosh2n-  - sin2£. 

r»  1 1 c 1 .1 


R.H.S.  > 


The  sign  and  the  magnitude  of  the  right-hand  side  of  equation 

(8)  is  used  to  determine  the  applicable  .,  0.  . relationship, 

equations  (5)  to  (7).  For  example,  if  using  the  most  recently 

computed  values  of  at  the  four  points  surrounding  i,  j yields  a 

R.H.S.  corresponding  to  equation  (18),  this  means  that  4/?  ^ >1.0 

and  the  new  value  of  . is  obtained  from: 

1 

R.H.S.  + (Ct+C.+Co+C,).  k„C  /k  C . (19) 


,j  (cosh2n.  - s i n 2 £ j ) / AFo  + ( (^ +C2+C3+C4)  k4Cps/k$CpJl 
The  new  value  of  ip^  . is  compared  with  the  old  guess  value 

■ * J 

to  see  if  it  satisifies  a convergence  criterion,  | . (new)  - 
, j (old ) I < 0.005. 

This  process  is  then  repeated  throughout  the  mesh  in  a definite 
order  a number  of  times  for  each  time  step  until  the  convergence 
criteria  is  met,  after  which  the  next  time  step,  ! is  calculated 
using  the  same  procedure. 
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V.  RESULTS  AND  DISCUSSION 


A.  Uniform  Heat  Flux 


The  analytical  solution  to  the  thermal  fields  in  a semi- 
infinite solid  substrate  subjected  to  a uniform  heat  flux  over  a 
circular  region  on  its  bounding  surface  is  given  by  Carslaw  and 
Jaeger  (5).  The  solution  may  be  simplified  to  give  the  temperature 
at  the  center  of  the  circular  region,  T(0,0).  One  form  of  this 
expression  is: 


k$  [T( 0 ,0)  - Tq] 


2/aTt 


qa 


/iT  a 


[1  - /rT  ierfc 


2/^t 


] 


(20) 


Equation  (20)  was  used  to  calculate  the  curve  associated  with 
the  vertical  axis  on  the  right  side  of  Figure  2.  This  plot  shows 
that  for  a given  radius  of  the  circular  region,  a,  there  is  a mini- 
mum heat  flux  required  if  the  center  of  the  circular  region  is  to 
reach  a given  temperature , e . g . the  melting  point  of  the  substrate.  That 


is,  for  small  values  of  a/ 2/a  t,  large  Fourier  numbers  at/a,  the 

S 5 


left  hand  side  of  equation  (20)  approaches  one  - the  temperature, 
T(0,0),  approaches  its  maximum  steady  state  value.  For  example,  in 


aluminum*  the  product  of  absorbed  heat  flux  and  the  radius  of  the 

.5 


circular  region  qa  ^ 1.45  x 10  W/m  if  the  center  of  the  circular 
region  is  to  ever  reach  the  melting  point. 

Similarly  one  may  anticipate  that  specific  qa^  values  could  be 
deduced  from  numerical  computations  for  the  initiation  of  surface 


^Properties  of  aluminum  used  in  the  calculations  are  listed  in 
Table  I. 


i 
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melting  and  for  any  surface  temperature  at  the  center  of  the  liquid 
zone  up  to  the  vaporization  temperature.  Figure  3 shows  a plot  of 
the  data  obtained  from  approximately  15  computer  runs  (over  two 
hundred  points).  Note  the  similarity  between  the  curve  associated 
with  the  vertical  axis  on  the  right  side  of  this  Figure  and  that 
in  Figure  2.  The  numerator  on  the  vertical  axis  includes  the  sum 
of  k£(T-TM),  ks(TM"To^  and  ksAHs«/CpJT  Tlie  latter  1S  used  here 
to  denote  the  equivalent  temperature  change  for  the  melting  of 
the  substrate.  The  curve  in  Figure  3 approaches  ^97%  of  the  steady 
state  value  for  small  a/2/a~t  values  which  is  sufficiently  accurate 
considering  other  approximations  used  in  the  computer  heat  flow 

5 

model.  For  an  aluminum  substrate,  minimum  values  of  qa  r 2.3  x 10 

5 

W/m  and  qa  ~ 4.2  x 10  W/m  are  deduced  from  this  curve  for  the 

initiation  of  surface  melting  and  for  a surface  temperature  T(0,0) 

5 

= Ty.  For  example  when  qa  = 2.86  x 10  the  calculated  maximum 
steady  state  temperature  at  the  center  of  the  molten  zone  is  1607  K. 
Equation  (20)  can  also  be  written  in  the  following  form: 


^ [T  (0,0)  - Tq] 
2 q/t 


,„  = !-/?  i erf c — - — 

* 2/^ 


! 


The  left  hand  side  of  equation  (21)  is  now  the  ratio  of  the 
surface  temperature  of  a semi -inf i ni te  solid  subjected  to  a uni- 
form heat  flux  over  its  entire  surface  to  the  temperature  T(0,0) 
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respectively.  In  either  case  the  curve  associated  with  the  verti- 
cal axis  on  the  left  side  of  Figure  2 shows  that  for  values  of 
a/2/a$t  > 2.0,  a$t/a  <:  0.0625,  heat  flow  is  essentially  unidirec- 
tional. A similar  plot  obtained  from  the  computer  heat  flow  model 
during  surface  melting  is  shown  in  Figure  3.  The  general  form  of 
the  two  curves  in  Figure  3 permit  the  deduction  of  the  following 
information  for  any  absorbed  heat  flux  and  radius  of  the  circular 
region.  First,  for  a given  value  of  £a  and  T(0,0)  one  may  estab- 
lish whether  lateral  heat  flow  in  the  substrate  can  be  ignored. 

As  example,  for  an  aluminum  substrate,  when  qa  ^ 1.4  x 10^  W/m 
heat  flow  is  essentially  in  one  dimension  for  T(0,0)  temperatures 
up  to  Ty.  Second,  the  tine  for  the  center  of  the  circular  region  to 

reach  any  temperature  above  the  melting  point  of  the  substrate  can 

g 

be  deduced.  For  example,  for  an  absorbed  uniform  heat  flux  of  10 

2 

W/m  and  a spot  radius  a = 380pm,  T(0,0)  will  reach  1700K  in  2.1  x 
_ 3 

10  sec.  Furthermore,  the  maximum  temperature  achieved  at  the 
center  of  the  circular  region  would  be  2238K. 

Figure  4 shows  the  general  relationship  developed  from  num- 
erical computations  between  temperature  at  the  center  of  the  cir- 
cular region  and  the  product  of  absorbed  heat  flux  and  melt  depth, 
qz,  for  an  aluminum  substrate,  z in  this  figure  denotes  melt  depth 
along  the  z-axis  at  the  center  of  the  circular  region  r = 0.  For 

5 

any  value  of  £ and  a,  for  which  the  product  qa  >,  2.3  x 10  W/m, 

time  for  T(0,0)  to  reach  a given  temperature  can  be  deduced  from 

Figure  3.  The  corresponding  melt  depth  along  the  z-axis  can  then 

be  found  from  Figure  4.  For  example,  for  uniform  heat  fluxes  of 
8 2 9 2 

3 x 10  W/m  and  10  W/ra  absorbed  over  circular  regions  with 
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radii  of  = 1587ym  and  s 476ym,  respectively  (qa  - 4.76  x 10  W/m) 
the  corresponding  times  and  melt  depths  during  melting  of  an  alum- 

_ O 

inum  substrate  when  T(0,0)  = are  t = 5.4  x 10  sec  and  t = 4.86 
x 1CT  sec  from  Figure  3 and  z = 933ym  and  z = 280ym  from  Figure  4, 
respectively.  Similar  information  can  also  be  obtained  from  these 
two  figures  for  any  T(0,0)  temperature  between  and  Tv.  While 
the  data  plotted  in  Figure  3 can  be  used  for  any  substrate  material, 
the  curve  in  Figure  4 is  specifically  calculated  for  an  aluminum 
substrate.  Similar  curves  can  readily  be  generated  for  other  sub- 
strate materials  using  the  equations  and  computations  procedure 
outlined  in  the  previous  sections.  As  noted  above,  for  values  of 
qa  > 1.4  x 10®  W/m  heat  flow  in  the  aluminum  substrate  is  essenti- 
ally one-dimensional  and  the  corresponding  times  and  melt  depths 
from  Figures  3 and  4 would  be  identical  to  those  previously  de- 
scribed (l)when  the  physical  and  thermal  properties  of  the  substrate 
are  the  same. 

An  important  finding  from  this  work  is  that  for  a given  value 
of  T ( 0 , 0 ) at  the  center  of  the  circular  region  the  dimensionless 
temperature  distribution  in  the  substrate  material  during  melting 
and  solidification  is  identical  for  any  combination  q^  and  a^  as 
long  as  qa  = constant.  Figure  5 shows  the  shape  and  location  of 
several  isotherms,  including  the  liquid-solid  interface,  in  an 
aluminum  substrate  material  subjected  to  a uniform  heat  flux  q 
over  a circular  region  of  radius  a^  where  the  product  qa  = 4.76  x 
105  W/m.  Note  that  these  are  the  isotherms  at  the  instant  T(0,0) 
reaches  the  vaporization  temperature  of  aluminum  Ty  = 2723K  and  the 
axis  used  are  made  dimensionless  dividing  them  by  the  radius  of 
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the  circular  region.  Identical  dimensionless  plots  were  obtained 
from  several  computer  runs  where  £ and  £ values  were  varied  over 
large  ranges  while  the  product  £a  was  kept  constant. 

Figure  6 shows  the  dimensionless  plots  of  the  location  and 

5 

shape  of  the  liquid-solid  interface  for  qa  = 4.76  x 10  W/m  during 
melting  and  subsequent  solidification  of  an  aluminum  substrate. 

In  these  simulations  the  heat  flux  was  removed,  PAV  in  equation  (15) 
was  set  equal  to  zero,  as  soon  as  T ( 0 , 0 ) reached  the  vaporization 
temperature.  As  in  the  previous  one-dimensional  heat  flow  model 
(1)  melting  of  the  substrate  continues  until  heat  flux  in  the 
liquid  and  the  solid  at  the  liquid-solid  interface  become  equal  - 
the  interface  velocity  approaches  zero.  Decreasing  values  of 
a/2/a~T,  increasing  time,  denote  the  successive  positions  of  the 
liquid-solid  interface  during  the  melting  and  the  subsequent  sol- 
idification sequence.  The  curve  in  Figure  6 identified  by 
a/2/astmax  shows  the  shape  and  position  of  the  liquid-solid  inter- 
face when  the  melt  along  the  z-axis  reaches  its  maximum  depth.  It 
is  interesting  to  note  that  due  to  the  lower  temperature  gradients 
in  the  liquid  near  the  perphery  of  the  molten  region  solidification 
in  this  region  commenced  while  melting  in  the  central  region  was 
still  continuing.  The  higher  lateral  heat  loss  thus  results  in 
liquid-solid  interface  shapes  closer  to  that  of  a spherical  geometry 
during  the  solidification  sequence.  Finally,  both  dimensionless  time 
a/2/aHt,  and  liquid-solid  interface  shapes  and  velocities  remained 
the  same  in  these  types  of  dimensionless  plots  as  long  as  the  product 
qa  was  kept  constant. 

Figures  7 and  8 show  the  effect  of  increasing  the  product  £a_  on 
the  geometry  and  the  location  of  the  liquid-solid  interface  at  the 
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instant  the  temperature  T(0,0)  reaches  2100K  and  2723K,  respectively. 
The  curves  for  qa  - 1.9  x 10^  W/m  are  in  the  range  where  heat  flow 
can  be  assumed  to  be  essentially  one-dimensional  and  the  shape  of  the 
isotherms  verify  this  fact.  On  the  other  hand,  with  decreasing  val- 
ues of  £a,  lateral  heat  flow  becomes  significant  and  the  liquid- 
solid  interfaces  assume  shapes  closer  to  the  oblate  spheroidal  coor- 
dinate illustrated  in  Figure  1. 

Liquid-solid  interface  velocities  during  melting  and  sub- 
sequent solidification  of  an  aluminum  substrate  along  the  z-axis 
and  the  r-axis  as  a function  of  £a  are  shown  in  Figures  9 and  10. 
These  data  show  that  for  a given  value  of  £_a  the  interface  velocity, 
R,  is  directly  proportional  to  the  absorbed  heat  flux  during  both 
melting  and  subsequent  solidification.  Furthermore,  for  any  given 
value  of  £a  interface  velocities  along  the  r-axis,  on  the  specimen 
surface,  are  higher  than  those  along  the  z-axis. 

Figures  9 and  10  also  show  that  in  the  one-dimensional  heat 

flow  range,  qa  » 1.4  x 10^  W/m,  liquid-solid  interface  velocities 

are  essentially  identical  to  that  previously  reported  (1).  Inter- 

★ 

face  velocities  rapidly  increase  to  a maximum  value  once  surface 
melting  commences  and  drop  off  to  zero  sometime  after  T(0,0)  = Ty 
and  the  heat  flux  is  removed.  During  subsequent  solidification 
interface  velocities  increase  with  distance  solidified  but  always 
remain  lower  than  that  during  melting.  For  example,  for  an  absorbed 
heat  flux  of  q = 10^  W/m^  and  qa  = 1.905  x 10^  W/m  (a  = I905ym) 

*the  rapid  increase  after  the  initiation  of  melting  occurs  over  very 
small  fractional  distances  and  is  not  shown  in  these  Figures. 
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the  calculated  interface  velocities  along  the  z-axis  midway  during 
the  melting  and  the  solidification  sequences  (z/zmax  = 0.5)  are 
^0 .17  m/sec  and  ^0.15  m/sec,  respectively.  The  corresponding  inter- 
face velocities  in  the  radial  direction,  on  the  specimen  surface, 
at  r/rm,v  = 0*5  are  °°  3m/sec  and  n,  0.32m/sec,  respectively.  In- 
creasing  the  absorbed  heat  flux  by  a factor  of  five  while 
the  product  qa^  is  kept  constant  (a  = 381  pm)  would  result  in  a 

corresponding  factor  of  five  increase  in  the  interface  vel- 
ocities calculated  above. 

5 6 

For  values  of  2.3  x 10  < qa  < 1.4  x 10  W/m  the  liquid-solid 

interface  velocities  along  the  z-direction  are  always  lower  during 
melting  than  during  solidification.  This  is  simply  from  the  in- 
creasing lateral  heat  loss  due  to  the  increased  divergency  of  heat 
flow  with  decreasing  values  of  qa  in  this  range.  However,  for 

5 

lower  values  of  qq,  e.g.  qa  = 4.76  x 10  W/m,  interface  velocities 
along  the  r-axis  are  approximately  the  same  during  melting  and 

g 

solidification.  For  example,  at  an  absorbed  heat  flux  of  q = 10 
W/m2  and  qa  = 4.76  x 105W/m(a  = 476pm)  the  calculated  interface  vel- 
ocities along  the  z-axis  midway  during  the  melting  and  the  solid- 
ification sequences  are  ^0.063m/s  and  ^0.24m/s,  respectively.  The 
corresponding  values  in  the  radial  direction,  from  Figure  10,  are 
^0.68m/s  and  ^0.57m/s,  respectively. 

No  attempt  is  made  here  to  quantitatively  describe  the  more 
complex  relationships  between  absorbed  heat  flux  and  the  ratio  of 
temperature  gradient  in  the  liquid  at  the  liquid-solid  interface 
and  growth  velocity,  G^/R . However,  qualitative  relationships  are 
easily  deduced  from  Figures  4 to  10.  In  general,  for  a given  value 
of  qq  temperature  gradient  at  the  liquid-solid  interface  increases 
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proportionally  with  increasing  values  of  absorbed  heat  flux,  £. 

On  the  other  hand,  since  interface  velocities  also  increase 
with  increasing  values  of  q,  the  ratio  G^/R  is  not  affected.  The 
maximum  values  of  GL/R  during  solidification  are  noted  at  the  ini- 
tial stages  of  solidification,  when  R is  almost  zero  and  the  temp- 
erature gradient  in  the  liquid  is  high.  Furthermore,  superheat 
in  the  liquid  is  rapidly  consumed  while  interface  velocities  in- 
crease with  increasing  distance  solidified  - resulting  in  a rapid 
decrease  in  GL/R.  Figures  6,  9 and  10  also  show  that  temperature 
gradients  in  the  liquid  are  higher  and  interface  velocities  are 
lower  along  the  z-axis  during  solidification,  therefore,  the  ratio 
G^/R  assumes  its  maximum  value  along  this  axis. 


B . Gaussian  Absorbed  Heat  Flux 

In  order  to  systemmatical ly  analyze  the  effect  of  Gaussian 
absorbed  heat  flux  distributions  on  the  surface  melting  and  solid- 
ification parameters  and  relate  same  to  the  previous  case  of  uni- 
form absorbed  heat  flux,  the  relationship  between  the  two  heat 
flux  distributions  shown  in  Figure  11  was  first  developed.  Inside 
the  circular  region  the  oblate  spheroidal  coordinate  n = 0,  Figure 
1.  Equation  (16)  for  case  (b)  can  thus  be  written  in  its  equiv- 
alent form  in  cylindrical  coordinates: 


% e 


■2r2/a2 


(22) 


If  we  assume  that  the  total  power  absorbed,  Q,  in  the  circu- 
lar region  is  identical  for  the  uniform  and  Gaussian  heat  flux 
distributions,  as  shown  in  Figure  11,  then  the  following  relation- 
ship is  deduced: 
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^uniform  ~ 2.313 

where  qQ  i s the  absorbed  heat  flux  at  the  center  of  the  circular 
region  and  will  be  referred  to  hereafter  as  the  maximum  absorbed 
heat  flux. 

The  equations  and  solution  methods  described  in  previous 
sections  were  used  to  generate  data  for  the  rapii  melting  and 
subsequent  solidification  of  an  aluminum  substrate  material. 

The  relationship  between  melt  depth  along  the  z-axis,  maximum 
absorbed  heat  flux,  qQ,  and  temperature  at  the  center  of  the  cir- 
cular region  is  shown  in  Figure  4.  The  data  indicate  that  for  a 
given  temperature  T(0,0)  and  total  power  absorbed,  Q,  melt  depths 
in  the  Gaussian  heat  flux  distribution  are  ^2.313  times  shallower 
than  that  obtained  in  the  case  of  the  uniform  absorbed  heat  flux. 

Computer  calculations  also  showed  that  dimensionless  temp- 
erature distributions  during  surface  melting  and  subsequent  solid- 
ification are  identical  if  the  product  qQa  is  kept  constant  while 
qQ  and  a are  varied.  Figure  12  shows  the  shape  and  location  of 
several  isotherms,  including  the  liquid-solid  interface  (T  * 933K), 
during  surface  melting  of  an  aluminum  substrate  subjected  to  a 
Gaussian  absorbed  heat  flux  distribution.  These  are  the  isotherm 
shapes  at  the  instant  when  T(0,0)  = 2723K. 

Figure  13  shows  the  effect  of  increasing  the  product  qQa/2.313 
on  the  geometry  and  the  location  of  the  liquid-solid  interface  at 
the  instant  when  T(0,0)  reaches  the  vaporization  temperature  of 
the  aluminum  substrate.  Increasing  the  product  qQa  results  in 


shallower  melt  depths  as  anticipated.  It  should  be  noted  that  the 


center  of  the  circular  region  reaches  the  vaporization  tempera- 
ture of  the  substrate  for  qQa/2.313  values  lower  than  the  criteria 
developed  for  the  uniform  absorbed  heat  flux,  qa  > 4.2  x 105  W/m. 
This  is  simply  due  to  the  higher  absorbed  heat  flux  at  the  center 
of  the  circular  region  in  the  Gaussian  distribution. 

Figures  14  and  15  compare  the  shapes  and  locations  of  the 
liquid-solid  interfaces  between  uniform  and  Gaussian  heat  flux 
distributions  when  T(0,0)  reaches  1420K  and  2723K,  respectively. 

As  indicated  above,  for  a given  total  absorbed  heat  flux  equation 
(23)  holds  and  the  center  of  the  circular  region  heats  up  much  more 
rapidly  when  the  heat  flux  has  a Gaussian  distribution.  The  melt- 
depth  along  the  z-axis  is  shallower,  while  the  melt  width  is  smaller 
due  to  the  diminishing  absorbed  heat  flux  with  increasing  r in  the 
circular  reaion. 

Figure  16  shows  dimensionless  plots  of  the  location  and  shape 

5 

of  the  liquid-solid  interface  for  qQa/2.313  = 4.76  x 10  W/m  during 
melting  and  subsequent  solidification  of  the  surface  of  an  aluminum 
substrate.  As  before,  the  heat  flux  was  removed  as  soon  as  T(0,0) 
reached  the  vaporization  temperature.  Note  the  similarity  between 
Figure  16  and  6,  the  total  absorbed  power  over  the  circular  region 
was  the  same  in  the  two  cases.  The  higher  maximum  heat  flux  in 
the  Gaussian  distribution  results  in  higher  temperature  gradients, 
shallower  melt  depths  and  shorter  times  - larger  values  of  a/2/a^T. 
For  example,  for  a total  absorbed  heat  flux  Q - 0.71  KW  over  a 
circular  region  of  radius  a - 476pm  the  times  for  the  interface  to 
reach  its  maximum  melt  depth  for  uniform  and  Gaussian  heat  flux 
distributions  are  5.9  x 10~3  second  6.25  x 10"^sec,  respectively. 


Liquid-solid  interface  velocities  during  melting  and  sub- 
sequent solidification  of  an  aluminum  substrate  along  the  z-axis 
and  the  r-axis  as  a function  of  qQa/2. 31 3 are  shown  in  Figures  17 
and  18,  respectively.  Note  that  in  these  two  figures  the  inter- 
face velocity,  R,  is  divided  by  qQ/2.313  - the  equivalent  heat 
flux  in  the  uniform  distribution.  These  data  again  show  that  for 
a given  value  of  qQa/2.313  the  interface  velocity  is  directly  pro- 
portional to  the  absorbed  heat  flux,  and  that  velocities  along  the 
r-axis  are  higher  than  those  along  the  z-axis.  Furthermore,  as 
the  product  of  qQa/2.313  decreases,  heat  flow  lines,  during  solid- 
ification, become  more  divergent  resulting  in  a corresponding  rel- 
ative increase  in  interface  velocities  along  the  z-axis.  Finally, 


for  a given  total  power  absorbed, the  calculated  interface  velocities 
along  the  z-axis  in  the  Gaussian  distribution  are  higher  than  those 
for  the  uniform  absorbed  heat  flux.  For  example,  at  an  absorbed 
heat  flux  of  qQ/2 . 313  = 109  W/m2  and  qQa/ 2 .31  3 = 4.76  x 1 05  W/m 
(a  = 476pm)  the  interface  velocities  from  Figure  17  during  melting 
and  solidification  are  0.27  m/s  and  0.47  m/s,  respectively.  The 
corresponding  values  for  the  same  total  power  absorbed  uniformly 
over  the  circular  region  were  previously  noted  to  be  ^0.063  m/s 
and  'v/  0.24  m/s,  respectively.  The  same  observation  holds  true  when 
interface  velocities  in  the  radial  direction  are  compared.  Figures 
18  and  10. 
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VI.  SUMMARY 

The  concepts  of  an  enthalpy  model  recently  proposed  are  used 
to  develop  the  heat  flow  equations  in  oblate  spheroidal  coordinates 
for  melting  and  solidification  of  the  surface  of  a semi -i nf i ni te 
solid  subjected  to  a variable  heat  flux  in  a circular  region  on 
its  bounding  surface.  While  numerical  solutions  are  presented  for 
step  function  uniform  and  Gaussian  heat  flux  distributions,  the 
equations  and  numerical  technique  developed  can  be  applied  to  heat 
fluxes  that  are  both  a function  of  position  and  time.  General 
trends  established  relate  the  important  melting  and  solidification 
parameters  to  absorbed  heat  flux,  radius  of  the  circular  region 
and  time  using  the  thermophysical  properties  of  an  aluminum  sub- 
strate material.  No  account  is  taken  of  the  changes  in  reflectiv- 
ity of  substrate  with  temperature  or  convection  in  the  molten 
region.  The  former  can  be  readily  included  in  the  generalized 
equations  when  experimentally  measured  values  become  available, 
while  the  latter  could  be  approximately  accounted  for  by  increasing 
the  thermal  conductivity  of  the  liquid. 


NOMENCLATURE 
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radius  of  the  circular  region,  m 


integration  constant 
specific  heat,  JKg’^K"^ 


fraction  liquid 


Fourier  number  (kct/pC  aj;  AFo 

b p 

dimensionless  time  step, 
specific  enthalpy,  JKg”^ 
heat  of  fusion,  JKg”^ 

thermal  conductivity,  Jm‘^s~^K"^ 

rate  of  heat  generation  per  unit  volume, 

Wm  "3 

_ ? 

absorbed  heat  flux,  Wm 
rate  of  total  absorbed  heat,  W 
radial  coordinate,  m 
time,  s 

time  for  the  surface  of  a material  to  reach  its 
melting  temperature, 

time  for  a material  to  achieve  a maximum  melt 
depth, 

time  for  surface  of  a material  to  reach  its  vapor- 
ization temperature, 

ambient  temperature  K 

melting  temperature  K 

surface  temperature  K 

vaporization  temperature  K 

volume,  m 

volume  of  discritized  domain 

axial  coordinate,  m 

thermal  diffusivity  (K/pCp),  m2  s"1 
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n oblate  spheroidal  coordinate 

0 dimensionless  temperature  variable,  equation  (4) 

C oblate  spheroidal  coordinate 

< p dimensionless  enthalpy  variable,  equation  (4) 

p density  Kg  m 


Subscripts 


1 * J 

z 

s 


nodal  point  subscripts  in  n and  £ directions, 
respecti vely 

liquid  region 

solid  region 


Superscri pts 

m time  level 


— - •- — - ..I., — • 
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LIST  OF  FIGURE  CAPTIONS 


Figure  1.  Two-dimensional  spheroidal  coordinate  system  and 
problem  geometry  for  the  rapid  melting  and  solid- 


ification of  a surface  layer. 

Figure  2.  Temperature  at  the  center  of  a circular  region  of 
a semi-infinite  solid  substrate  as  a function  of 
applied  uniform  absorbed  heat  flux,  radius  of  the 
circular  region  and  time.  The  vertical  axis  on  the 
left  gives  the  ratio  between  temperatures  in  one 
and  two-dimensional  heat  flow. 

Figure  3.  Temperature  at  the  center  of  the  liquid  zone  of  a 

semi-inf ini te  solid  substrate  during  surface  melting 
as  a function  of  uniform  absorbed  heat  flux,  radius 
of  the  circular  region  and  time.  The  vertical  axis 
on  the  left  gives  the  ratio  of  square  root  of  times 
in  one  and  two-dimensional  heat  flow. 

Figure  4.  Temperature  at  the  center  of  the  circular  region 

T(0,0)  over  which  a heat  flux  is  absorbed  in  an  alum- 
inum substrate  versus  the  product  of  the  absorbed  heat 
flux  and  melt  depth  in  the  axial  direction  away  from 
the  surface.  qQ  denotes  the  maximum  heat  flux  at  the 
center  of  the  Gaussian  heat  flux  distribution.  See 
Figure  11  and  equation  (23). 

Figure  5.  Shape  and  location  of  several  isotherms,  including  the 
liquid-solid  interface  (T  = 933K),  during  melting  of 
an  aluminum  substrate  subjected  to  a uniform  absorbed 
heat  flux  £ over  a circular  region  of  radius  a.  The 
isotherms  are  drawn  at  the  instant  when  T(0,0)=TV=2723K, 
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Figure  6.  Liquid-solid  interfaces  during  melting  and  subsequent 
solidification  of  an  aluminum  substrate  subjected  to 
a uniform  absorbed  heat  flux  £ over  a circular  region 
of  radius  a.  The  curve  marked  T(0,0)=Tv  shows  the 
location  of  the  interface  at  which  time  the  heat  flux 
was  removed.  Decreasing  values  of  a / 2 /a  t denote  the 
successive  positions  of  the  liquid-solid  interface 
during  the  melting  and  the  subsequent  solidification 
sequence . 

Figure  7.  Liquid-solid  interfaces  during  melting  of  an  aluminum 

substrate  subjected  to  a uniform  heat  flux  £ over  a cir- 
cular region  of  radius  a.  The  interfaces  for  each  con- 
stant ££  are  shown  at  the  instant  the  center  of  the  cir- 
cular region  reaches  the  temperature  T( 0 , 0) =21 00K. 

Figure  8.  Liquid-solid  interfaces  during  melting  of  an  aluminum 

substrate  subjected  to  a uniform  heat  flux  £ over  a cir- 
cular region  of  radius  £.  The  interfaces  for  each 
constant  £a  are  shown  at  the  instant  the  center  of  the 
circular  region  reaches  the  vaporization  temperature. 

Figure  9.  The  ratio  of  liquid-solid  interface  velocity  to  absorbed 
heat  flux,  R/q,  as  a function  of  fractional  distance, 
z/zmav,  during  melting  and  solidification  of  an  aluminum 
substrate  subjected  to  a uniform  absorbed  heat  flux,  £ 
over  a circular  region  of  radius  £. 

Figure  10.  The  ratio  of  liquid-solid  interface  velocity  to  absorbed 

heat  flux,  R/q,  as  a function  of  fractional  distance, 

r/r  , during  melting  and  solidification  of  an  aluminum 
mdx 

substrate  subjected  to  a uniform  absorbed  heat  flux,  £ 
over  a circular  region  of  radius  £. 
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Figure  11. 


Figure  12. 


Figure  13. 


Figure  14. 


Figure  15. 


Figure  16. 


Relationship  between  uniform  and  Gaussian  absorbed 
heat  fluxes  when  the  total  power  absorbed  over  the 
circular  region  is  identical. 

Shape  and  location  of  several  isotherms,  including 
the  liquid-solid  interface  (T=933K)  during  melting 
of  an  aluminum  substrate  subjected  to  a Gaussian 
absorbed  heat  flux  distribution  over  a circular  region 
of  radius  a.  The  isotherms  are  drawn  at  the  instant 
when  T(0,0)=Tv=2723K. 

Liquid-solid  interfaces  during  melting  of  an  aluminum 
substrate  subjected  to  a Gaussian  absorbed  heat  flux 
distribution  over  a circular  region  of  radius  a.  The 
interfaces  for  each  constant  qQa/2.313  are  shown  at  the 
instant  the  center  of  the  circular  region  reaches  the 
vaporization  temperature. 

Comparison  between  the  shapes  and  locations  of  the 
liquid-solid  interfaces  during  melting  of  an  aluminum 
substrate  subjected  to  uniform  and  Gaussian  absorbed 
heat  flux  distributions  at  the  instant  T(0,0)=1420K. 
Comparison  between  the  shapes  and  locations  of  the 
liquid-solid  interfaces  during  melting  of  an  aluminum 
substrate  subjected  to  uniform  and  Gaussian  absorbed 
heat  flux  distributions  at  the  instant  T(0,0)*Tv. 
Liquid-solid  interfaces  during  melting  and  subsequent 
solidification  of  an  aluminum  substrate  subjected  to 
a Gaussian  absorbed  heat  flux  distribution  over  a cir- 
cular region  of  radius  a^. 


Figure  17. 


The  ratio  of  liquid-solid  interface  velocity  to 
absorbed  heat  flux  as  a function  of  fractional 


distance,  z/z  , during  melting  and  solidifi- 
cation  of  an  aluminum  substrate  subjected  to  a 
Gaussian  absorbed  heat  flux  distribution  over  a 
circular  region  of  radius  a. 

Figure  18.  The  ratio  of  liquid-solid  interface  velocity  to 
absorbed  heat  flux  as  a function  of  fractional 
distance,  r/rmax»  during  melting  and  solidifi- 
cation of  an  aluminum  substrate  subjected  to  a 
Gaussian  absorbed  heat  flux  over  a circular 
region  of  radius  a. 
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TABLE  I 

SUMMARY  OF  PROPERTIES  AND  SOME  CALCULATED 
CONSTANTS  FOR  THE  ALUMINUM  SUBSTRATE 


C * = 1048 
ps 

Cpt*  = 1086 


AH 


Si 

* 


V 
. * 


= 228 


k**  = 108 

TM  = 933 
= 2723 


J Kg”1  K-1,  Specific  Heat  of  Solid 
J Kg”1  K"1,  Specific  Heat  of  Liquid 


= 3.95  x 103  J Kg 


-1 


Latent  Heat  of  Fusion 


Wm”1  K"1,  Thermal  Conductivity  of  Solid 
Wm"1  K"1,  Thermal  Conductivity  of  Liquid 
K,  Melting  Temperature 
K,  Vaporization  Temperature 


-5  2 -1 

= 8.55  x 10  m sec  , Thermal  Diffusivity  of  Solid 


»**  -5  2 -I 

i = 3.91  x 10  m sec  , Thermal  Diffusivity  of  Liquid 


p = p***  = p***  = 2545  Kg  m" 


1-D 


= 1.5  x 10' 


. -2 
J m sec 


Dens i ty 
■1/2 


9^ 


= 4.5  x 10^  J m”^  sec”1^ 


1-D 


q/t, 


max 


1-D 


= 5.5  x 10^  J m”^  sec-1/^ 


* Averaged  from  298K  to  Melting  Point 

**  Averaged  from  Melting  Point  to  Vaporization  Temperature 

***  Average  of  Liquid  and  Solid 
1-D  Refers  to  One-Dimensional  Heat  Flow 
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Ficure  2.  Temperature  at  the  center  of  a circular  region  of 
a semi -i nf ini te  solid  substrate  as  a function  of 
applied  uniform  absorbed  heat  flux,  radius  of  the 
circular  region  and  time.  The  vertical  axis  on  the 
left  gives  the  ratio  between  temperature  in  one  and 
two-dimensional  heat  flow. 
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Figure  3.  Temperature  at  the  center  of  the  liquid  zone  of  a 

semi-inf inite  solid  substrate  during  surface  melting 
as  a function  of  uniform  absorbed  heat  flux,  radius 
of  the  circular  region  and  time.  The  vertical  axis 
on  the  left  gives  the  ratio  of  square  root  of  times 
in  one  and  two-dimensional  heat  flow. 
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Figure  4. 


Temperature  at  the  center  of  the  circular  region 
T(0,0)  over  which  a heat  flux  is  absorbed  in  an  alum- 
inum substrate  versus  the  product  of  the  absorbed  heat 
flux  and  melt  depth  in  the  axial  direction  away  from 
the  surface.  q0  denotes  the  maximum  heat  flux  at  the 
center  of  the  Gaussian  heat  flux  distribution.  See 
Figure  11  and  equation  (23). 
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Figure  6.  Liquid-solid  interfaces  during  melting  and  subsequent 
solidification  of  an  aluminum  substrate  subjected  to 
a uniform  absorbed  heat  flux  q over  a circular  region 
of  radius  a.  The  curve  marked  T(0,0)=T  shows  the 
location  o7  the  interface  at  which  timevthe  heat  flux 
was  removed.  Decreasing  values  of  a/2/a.t  denote  the 
successive  positions  of  the  liquid-sol idsinterface 
during  the  melting  and  the  subsequent  solidification 
seouenco . 
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Figure  7.  Liquid-solid  interfaces  during  melting  of  an  aluminum 

substrate  subjected  to  a uniform  heat  flux  £ over  a cir 
cular  region  of  radius  a^.  The  interfaces  for  each 
constant  £a  are  shown  at  the  instant  the  center  of  the 
circular  region  reaches  the  temperature  T(0,0)=2100K. 


Figure  8. 


Liquid-solid  interfaces  during  melting  of  an  alum- 
inum substrate  subjected  to  a uniform  heat  flux  £ 
over  a circular  region  of  radius  a^.  The  interfaces 
for  each  constant  3a  are  shown  at  the  instant  the 
center  of  the  clrcuTar  region  reaches  the  vaporiza- 
tion temperature. 
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Figure  9.  The  ratio  of  liquid-solid  interface  velocity  to  absorbed 
heat  flux,  R/q,  as  a function  of  fractional  distance, 
z/2max»  during  melting  and  solidification  of  an  alum- 
inum substrate  subjected  to  a uniform  absorbed  heat 
flux,  ^ over  a circular  region  of  radius  a^ 
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The  ratio  of  liouid-solid  interface  velocity  to  absorbed 
heat  flux,  R/q,  as  a function  of  fractional  distance, 
r/rmax»  during  melting  and  solidification  of  an  alum- 
inum substrate  subjected  to  a uniform  absorbed  heat 
flux,  q over  a circular  region  of  radius  &. 


Figure  10. 
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Figure  12.  Shape  and  location  of  several  isotherms,  including 
the  liquid-solid  interface  (T=933K)  during  melting 
of  an  aluminum  substrate  subjected  to  a Gaussian 
absorbed  heat  flux  distribution  over  a'  circular  region 
of  radius  a.  The  isotherms  are  drawn  at  the  instant 
when  T(0,0?=TV=2723K. 
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Figure  13.  Liquid-solid  interfaces  during  melting  of  an  aluminum 
substrate  subjected  to  a Gaussian  absorbed  heat  flux 
distribution  over  a circular  region  of  radius  a.  The 
interfaces  for  each  constant  q0a/2.313  are  shown  at 
the  instant  the  center  of  the  circular  region  reaches 
vaporization  temperature. 
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liquid-solid  interfaces  during  melting  of  an  alum- 
inum substrate  subjected  to  uniform  and  Gaussian 
absorbed  heat  flux  distributions  at  the  instant 
T(0,0)=1420K. 


Figure  15.  Comparison  between  the  shapes  and  locations  of 
the  liquid-solid  interfaces  during  melting  of 
an  aluminum  substrate  subjected  to  uniform  and 
Gaussian  absorbed  heat  flux  distributions  at  the 
instant  T(0,0)=T  . 
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Figure  16.  Liquid-solid  interfaces  during  melting  and  subse- 
quent solidification  of  an  aluminum  substrate  sub- 
jected to  a Gaussian  absorbed  heat  flux  distribution 
over  a circular  region  of  radius  a. 


Figure  17.  The  raio  of  liquid-solid  interface  velocity  to 
absorbed  heat  flux  as  a function  of  fractional 
distance,  z/zmaX,  during  melting  and  solidifi- 
cation of  an  aluminum  substrate  subjected  to  a 
Gaussian  absorbed  heat  flux  distribution  over  a 
circular  region  of  radius  a. 
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Figure  18.  The  ratio  of  liquid-solid  interface  velocity  to 
absorbed  heat  flux  as  a function  of  fractional 
distance,  r/rmaX,  during  melting  and  solidifi- 
cation of  an  aluminum  substrate  subjected  to  a 
Gaussian  absorbed  heat  flux  over  a circular  region 
of  radius  a. 


